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Abstract 

2 Using the recently developed effective field theory of inflation, we argue that the size and the shape 
^ of the non-Gaussianities generated by single-field inflation are generically well described by two 
^ pa_: .^which ch^cterizes .he of th. .ignal tha. is peaKed on cuila.cal co„- 

figurations, and f^^ , which instead characterizes the size of the signal which is peaked both 

> 



on equilateral configurations and flat-triangle configurations (with opposite signs). The shape of 
non-Gaussianities associated with f'^]l^°^' is orthogonal to the one associated to /^^^'', and former 
^ analysis have been mostly blind to it. We perform the optimal analysis of the WMAP 5-year data 
If^ for both of these parameters. We find no evidence of non-Gaussianity, and we have the following 
O constraints: -125 < fj^f' < 435, -369 < f^*^°^- < 71 at 95% CL. We show that both of these 
constraints can be translated into limits on parameters of the Lagrangian of single-field inflation. 
^ For one of them, the speed of sound of the inflaton fluctuations, we find that it is either bounded 
to be Cg > 0.011 at 95% CL. or alternatively to be so small that the higher-derivative kinetic term 
^ dominate at horizon crossing. We are able to put similar constraints on the other operators of the 
inflaton Lagrangian. 
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1 Introduction 



The interest in measuring a possible non-Gaussianity of the primordial density perturbation 
has grown in the last few years due mainly to two reasons: on the one hand, new experiments 
have finally become sensitive enough to measure even relatively small deviations from Gaus- 
sian statistics; on the other hand, new theoretical insights have shown that while standard 
slow-roll inflation where the inflaton has a canonical kinetic term predicts an undetectable 
level of non-Gaussianity [1], many other inflationary models can give rise to a level of non- 
Gaussianity already detectable by current experiments. Furthermore it has been shown [2, 3] 
that a detection of non-Gaussianity would carry an unprecedented amount of information on 
the dynamics that drove inflation, allowing us to explore the interaction Lagrangian of the 
inflaton. 

The current constraints on the level of non-Gaussianity of the primordial density pertur- 
bations already imply that the deviation from a Gaussian spectrum is at most at the percent 
level. This means that the most natural observables to look at are ones which are zero in the 
presence of an exactly Gaussian spectrum. This is the case for the bispectrum 

i^A'^ks) = {'27^rS^'\J2^^)F{k„h,h) , (1) 

i 

where $ is the Newtonian potential. Because of translation invariance, the sum of the spatial 
/c's must form a closed triangle. In the case of the non-Gaussianities produced by inflation, 
approximate scale invariance implies that, under rcscaling of the momenta, the F above scales 
approximately as where k is the typical size of one of the momenta. Finally, rotation 

invariance forces F to be only a function of the shape of the triangle (which one can think as 
being parameterized by the ratios of the two lowest /c's with respect to the highest k in the 
triangle) . In some models F can deviate from scale invariance in a way that might have some 
measurable effects if one looks at non-Gaussianities on very different scales [4]. At present 
this effect is not very important for Cosmic Microwave Background (GMB) experiments given 
the lack of a detection. Therefore we will neglect this in what follows. 

So far, all analyses both in the CMB and in Large Scale Structure (LSS) have concentrated 
just on two kinds of non-Gaussianities. The first one is the so called local type, in which $ is 
given by [5] 

= + fj^t (<I>(^)(f)2 - , (2) 

where ^^^^ is a Gaussian stochastic variable. determines the size of the non-Gaussianities, 
and is the parameter being constrained. In momentum space this kind of non-Gaussianity 
implies that the signal is peaked on 'squeezed' triangles where one of the sides is much smaller 
than the other two. This kind of non-Gaussianity is predicted in some multi-field infiationary 
models, see for example [6, 7], or in the recently proposed new bouncing cosmology [8]. There 
is a theorem that states that the non-Gaussian signal in the squeezed limit that can be 
produced by single-field inflation, has to be proportional to the tilt of the power spectrum 
[1, 9, 10]. This means that, unless a strong deviation from scale invariance in the power 
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spectrum is found, detection of such a non-Gaussian signal would definitively rule out all 
single-field inflationary models. The current constraint on this kind of signal from the WMAP 
5yr data is given by [11]: 

-4 < /nI < 80 at 95% CL , (3) 

where the optimal estimator found in [12] was implemented with the pipeline developed in 
[13]. When combined with the constraint from LSS [14], this gives — 1 < < 63 at 95% 
CL. 

The second kind of non-Gaussianity that has been constrained is the so-called equilateral 
kind. In this case, the signal is peaked on equilateral-triangle configurations in Fourier space. 
This is the shape that is produced in some proposed single-field or multi-field infiationary 
models, such as DBI infiation [15, 16, 17], or Ghost infiation [18, 19], where the infiaton has 
large derivative interactions, or Trapped inflation [21], where production of massive particles 
slows down the infiaton. This shape is quite different from the local shape, and requires 
a dedicated analysis. The current best limits are given by the WMAP team, who applied 
the technique developed in [12] to the WMAP Syr data: -151 < f^f' < 253 at 95% CL 
[22]. In this paper we will improve the above limit by performing the optimal analysis for 
the same parameter /^^''', using the estimator first presented in [12], and implementing it 
with the improved pipeline of [13]. Furthermore, this allows us also to marginalize over the 
amplitude of the templates of the galactic foregrounds, performing an improved treatment of 
the foreground with respect to the cleaning procedure applied in [22]. After implementing 
the optimal analysis, we find no evidence of non-Gaussianity: 

-125 < fPf- < 435 at 95% CL , (4) 

This constraint is obtained by assuming that other forms of primordial non-Gaussianities are 
absent. Though this limit is optimal, it is a little bit milder than the one claimed by the 
WMAP collaboration on the same set of data. As we explain in detail in sec. 4.2, we believe 
that previous analyses have underestimated the error ^• 

Improving the limit on /^^^^' is not however the main purpose of the present paper. Thanks 
to the development of the effective field theory of infiation [3] , we will be able to argue that the 
parameter space of the non-Gaussianities produced by the most general single-field models, 
where the infiaton fiuctuations have an approximate shift symmetry, is in reality larger than 
the one characterized by /^^'^'. It consists of any linear combination of two independent 
shapes: the equilateral one, and a new one that we call orthogonal. This orthogonal shape 
is peaked both on equilateral-triangle configurations and on flattened-triangle configurations 
(where the two lowest-Zc sides are equal exactly to half of the highest-/^ side). The sign in this 
two limits is opposite. The new shape is orthogonal to the equilateral one (we define a scalar 
product later in the text). We define f^^j^"^' to parameterize the amplitude of the orthogonal 
template, and we then argue that former analysis of the bispectrum for the local and the 

^Wc acknowledge Eiichiro Komatsu for collaboration regarding some comparisons of our code with the one 
used by the WMAP team. 
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equilateral shape have been quite insensitive to the orthogonal shape which therefore could 
be rather large in the data and still not have been detected. After performing the optimal 
analysis, we find no evidence of a non-zero fp^i^°^', and we obtain the constraint: 

-369 < < 71 at 95% CL . (5) 

In addition, thanks again to the effective field theory formalism, we are able to map the 
constraint we find on f^l'^' and f^]i^°^' into constraints on the coefficients of some operators of 
the Lagrangian for the inflationary fluctuations (under of course the assumption that there is 
only one relevant light degree of freedom during the inflationary phase). We believe that being 
able to put constraints or potentially measure the coefficients of the interaction Lagrangian in 
inflation is quite a remarkable fact, given that the energy scale at which this can happen could 
be huge. For the readers familiar with particle physics and collider physics, loosely speaking 
the limits we will put on the inflaton Lagrangian can be seen as the analogue of the constraints 
on the higher dimension operators that come from precision electroweak tests [23, 24, 25]. The 
only difference is that while in collider physics we might hope to produce some of the particles 
mediating new processes on shell at particle accelerators, it is very unlikely that we will be 
ever able to produce directly the inflaton particle. We believe this makes the formalism we are 
developing even more important. This comparison with collider physics is somewhat natural, 
as having non-Gaussianities in the sky, with their huge amount of information associated 
to the different triangular configurations, is very similar to having in an accelerator a cross 
section as a function of the angle of the outgoing particles. The inclusion of the orthogonal 
shape in the analysis allows us to put constraints on the parameter space of the interactions 
of the inflationary fluctuations. Since the presence of a particular kind of non-Gaussianity is 
associated to a small speed of sound of the fluctuations, we can constrain Cg to be larger 
than: 

> 0.011 at 95% GL . (6) 

or to be smaller than 

cs < 10-2 {d2 + d.f" , (7) 

In this last case the higher-derivative kinetic term, whose size is proportional to (^2 + dz), is 
important at horizon crossing and the non-Gaussianities depend on other coefficients. As we 
argue in the text, without the inclusion of the orthogonal shape into the analysis, it would 
have been impossible to set this bound without rather strong assumptions. We are then able 
to derive similar constraints on other three parameters of the inflaton Lagrangian. In some 
region of the parameter space, consistent inflationary models have a negative squared speed 
of sound for the fluctuations at horizon crossing. This leads to an exponential growth of 
the perturbations before horizon crossing. The analysis of the WMAP data shows that these 
models are practically ruled out at 95% GL. 

The paper is organized as follows. In sec. 2 we briefly review the effective held theory of 
inflation, and we show that there are two independent shapes that need to be analyzed. In 
sec. 3 we show how we can analyze the two shapes in the GMB data. In sec. 4, we actually 
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give the results of the analysis of the WMAP data and the constraints on the Lagrangian for 
the inflationary fluctuations. In sec. 5, we summarize our results. 



2 Effective Field Theory of Inflation and Shape of non- 
Gaussianities 

An effective field theory description for the fluctuations of the inflaton has been developed 
in [3]. By unifying in one description all single- field models, it allows us to explore in full 
generality the signature space of single-field inflation, in its broadest sense. In inflation there 
is a physical clock that controls when inflation ends. This means that time translations are 
spontaneously broken, and that therefore there is a Goldstone boson associated with this 
symmetry breaking. As usual, the Lagrangian of this Goldstone boson is highly constrained 
by the symmetries of the problem, in this case the fact that the spacctime is approximately 
de Sitter space, in the sense that <C 1. The Goldstone boson, that we can call tt, 

can be thought of as being equivalent, in standard models of inflation driven by a scalar field, 
to the perturbations in the scalar field 50. The relation valid at linear order is tt = 50/0, 
where is the speed of the background solution. We stress that the description in terms 
of the Goldstone boson tt is more general than this and it does not assume the presence 
of a fundamental scalar field. Although here and in the rest of the paper we often refer 
to the Lagrangian for the Goldstone boson as the Lagrangian for single-field inflation, this 
should be meant in the broadest sense that there is only one light relevant degree of freedom 
during the inflationary phase. It does not mean that the background solution is generated 
by a fundamental scalar field: the Lagrangian for the inflationary fluctuations expressed in 
terms of the Goldstone boson is universal and independent of the details through which the 
background solution is generated. 

The Goldstone boson tt is related to the standard curvature perturbation ( by the relation 
C = —Hit, which is valid at linear order and at leading order in the generalized slow roll 
parameters. The most general Lagrangian for the Goldstone boson is given by [3, 10]: 



HMl,.^ 2N-1/^ N2 HMl,.^ 2N A 253. 



+ 

where here we have assumed that the Goldstone boson is protected by an approximate shift 
symmetry that allows us to neglect terms where n appears without a derivative acting on 
it. Let us explain the symbols that appear in the above action. For the large level of non- 
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Gaussianities we will be interested in we can neglect the metric perturbations (see [3, 10] for 
a discussion about this approximation) and we can take the metric g^^j, to be unperturbed 
FRW: 

ds'^ = -dt'^ + a{tfdx^ . (9) 

H is the standard Hubble parameter, Cg is the speed of sound of the fluctuations, which, 
since the background is not Lorentz invariant, does not need to be cqiial to one. M is a free 
parameter with dimension of mass, while rfi, ^2, c^a and £3 are dimensionless parameters that 
as we will explain later are expected to be of order one ^. The dots represent higher derivative 
terms that in general give a negligible contribution to the observables, or quartic or higher 
order terms that we neglect because here we are interested in the 3-point function. 

The third line of (8) contains one time kinetic term and two spatial kinetic terms: the first 
two are standard ones whose coefficient is proportional to and the other one is a higher 
derivative term. There is also a higher-derivative cubic term. These terms become important 
when inflation happens very close to de Sitter space, in the sense that [3] 

Ml^\H\ < Max( diHM^/A, (^2 + d^)M^H^/{2 d;) ) . (10) 

This is so because the coefficient of (Sjvr)^ in the first line (and in some limits also the one 
of 7T^) is forced by the symmetries of the problem to be proportional to HMp^. As typical 
in the effective field theory approach, this fact appears so directly in the Lagrangian because 
the symmetries are fully exploited in its construction. In the limit in which the system is 
very close to de Sitter, this term can become subleading with respect to the ones of the third 
line of eq. (8), which are either a higher-derivative term, or a standard kinetic term, whose 
coefficient is however proportional to H. This happens when (10) is satisfied, in which case 
the system approaches Ghost Inflation [3, 10, 18, 19] or its generalization [3, 10, 20]. We will 
therefore neglect these three terms for the moment and we will come back to them later ^. 

We notice that, away from the near-de-Sitter limit, the interaction Lagrangian for the 
Goldstone boson contains, at leading order in derivatives, just two interaction operators: 
7r(c}j7r)^ and tt^. Notice that their coefficients are not fixed by the constraint of being close 
to a de Sitter phase; they can be large, and therefore induce large and detectable non- 
Gaussianities. These coefficients however are not completely free. The coefficient of 7r((9i7r)^ 
is fixed to be correlated with the speed of sound of the fluctuations (again, this appears so 
evidently, thanks to the effective held theory approach) . The coefficient becomes large in the 
limit of small speed of sound, which explains why a small speed of sound is associated to 
large non-Gaussianities. The coefficient of tt^ is instead dependent also on C3. We expect 
this number to be order one. This is the case for example in the UV complete model of DBI 
Infiation [15], where £3 = 3(1 — c^)/2 and Cg <^ 1. From the effective field theory point of 
view, when £3 is of order one, the strong coupling scale induced by the operator tt^ is the 

^Notice that £3 = c^Cs, where C3 is the parameter defined in [10]. The notation is chosen to be as close as 
possible to the one of [10]. 

^Strictly speaking, Cg is the speed of sound of the fluctuations only in the limit in which we neglect the 
term proportional to di, which is a good approximation when (10) is not satisfied. 
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same as the one induced by the operator 7r(9j7r)^ . The unitarity cutoff induced by 7r(5i7r)^ is 
in fact given by: 



AW-167r^^Pil^l(i:?^' (11) 
while the one due to tt^ is given by 

which are indeed of the same order for C3 of order one. It is also easy to estimate that if £3 
is order one, loop corrections renormalize the coefficients of the two operators only at order 
one level. This means that £3 of order one, and a small speed of sound c^, are technically 
natural ^. 

The three point function of the Newtonian potential $ has the usual form 

i^k.^A) = {2^f5^'\Y.h)F{hMM) . (13) 

i 

Here 

^ = 5C, (14) 

where C, is the curvature perturbation of comoving slices. This relationship is valid, at first 
order, out of the horizon, during matter domination. The 5— function comes from transla- 
tion invariance and it tells us that the 3-point function is a function of closed triangles in 
momentum space. For single-field infiation, F can be read off from: 

{^k,^hW = - (J) = (27r)^5(^HE i^HoMki, h, h) + F,s{h, h, h)) 

^ ' i 

(15) 

where F^(^q.^)2 is the shape generated by the operator 7r(9j7r)^, and F^s is instead the one 
generated by the operator tt^. In the limit in which we are far enough from de Sitter (in the 
sense of the inequality (10)), and in which we consider an approximate shift symmetry for 
the Goldstone boson, the resulting form of the non-Gaussianity is given by [10] ^: 

Fi,^Q^^y{h,k2,h) = (1 - ^) • ^1 (16) 
(24X36 - 8X2^X3=^X1 - 8X2^X12 + 22X3=^X1=^ - 6X2^X1^ + 2X1^) 

X3^X^^ ' 

^ /, , , X 20 / 1 \ , 3 2x .2 1 



3 V cjj ' 2 * X33Xi= 



^This discussion has assumed that the operator proportional to (^2 + ds) is ncghgible up the cutoff. This 
is the case only if (^2 + ds) < Cg. When this inequality is violated, (0^2 + ^3) niust be positive and, as we will 
later explain more in detail, the natural value of 63 gets scaled down by a factor of order cj^ jidi -\- ^3)^/^, 
which is clearly a negligible correction for [di -\- dz) of order one and for the values of Cg that are currently 
allowed by the data. 

^In a different and somewhat less general formalism, this expression was obtained also in [26], where it 
was already pointed out that the models considered in [26] admitted two independent shapes for the non- 
Gaussianities. 
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Here we have used that 



= (27r)3<5(3)(^i + 4)^ , (17) 



with 

9 



and e = —H/H"^. We have also defined 



^*- 2547^Jl^' ^^^^ 



= fci + k2 + A;3 , (19) 
K2 = (A;iA;2 + A:2A;3 + A;3A;i)'/' , 

= {klk2k^f'^ . 



We can use the standard definition of f^L 

k^ 



F{K k, k)^UL-—^, (20) 



to define 



= 1^ fl - ^ I . (21) 



324 V c2 
.3 10 / 1\ 3 



s 



and to write 



fcs) = /Jf'^^'Al (22) 



X 



17 

(24/^36 - 8K2'^Ks^Ki - 8K2'^Ki^ + 22K^^Ki^ - e/sTs^i^i^ + 2^1^) 



F^3(A;i,A;2,A;3) = 162/*>| ^ 



^337^13 

Notice that, quite remarkably, if wc wish to have a speed of sound smaller than one (a 
necessary condition for the existence of a Lorcntz invariant UV completion [27]), wc need to 
have fj,f^^' < 0. Even more importantly, the size and the shape of the three-point function is 
controlled by two parameters: the speed of sound Cg and the parameter 63, and the resulting 
non-Gaussian signal is a linear combination of two independent shapes. This means that 
if we fix the size of the overall non-Gaussianity there is a one-parameter family of shapes 
associated to the same amount of signal. In particular, for 63 of 0{1), a necessary condition 
for the 3-point function to be large is that the speed of sound is small ^. Given the current 
experimental sensitivity, we can therefore concentrate on this limit. As we let £3 vary keeping 
Cg fixed the size and the shape of the produced non-Gaussianity change. Concentrating on the 



^Another possibility is to be very close to the de Sitter limit. We will comment on this later. 
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shape, for very large or very small values of C3 the non-Gaussianity is dominated by one of the 
two operators, which both have a shape quite close to the so called equilateral kind. In that 
case the signal is concentrated on equilateral triangular configurations of the sort shown in 
the top-left panel of Fig. 1. However the two shapes are not identical which implies that there 
is a region at intermediate values of £3 where the shape of the resulting three-point function is 
completely different, being peaked on flat triangles where the size of the two lowest momenta 
is precisely half of the highest one. This can be seen in Fig. 1 where we show the shape of the 
non gaussianity for several values of the parameter £3. More quantitatively, in Fig. 2 we show 
the cosine between the generic shape F generated by single-field inflation (far from the near- 
de-Sitter limit) with F^(a.7r)2 as we let C3 vary. The cosine between two shapes was deflned 
in [2], and it represents a quantitative measure of the similarity and the correlation of the 
signals. Given two shapes F(i)(A;i, k2, k^) and F(2)(/ci, k2, ks), one flrst deflnes a 3-dimensional 
scalar product between the shapes as: 

^(1) ■ ^(2) = Yl h)Fi2){ku /C2, ks)/ {Pk.PkM , (23) 

j^physical 

where P(k) represents the power spectrum and kf^^^^'^^^ means that only the ^'s that form a 
triangle are included. The 3D cosine between two shapes is then deflned as 

Prom Fig. 2, we see that the cosine is very small around 63 ~ —5 for a region of approximately 
~ 10 — 20% of the parameter space. Here, somewhat arbitrarily, we consider that the natural 
parameter range for 63 is between -10 and 10, and that the region where the cosine is small is 
deflned as the region where this is smaller than 0.7. We realize that this is approximately the 
relevant number once we plot the scalar product of the shape with the local shape produced 
by multifleld models [6, 7] and by the new ekpyrotic universe [8]. This shape is plotted 
in Fig. 3. Although the shapes are clearly very different, we see that the cosine with the 
equilateral shape (£3 ~ 0) is approximately 0.4. In summary for roughly 10-20% of the 
natural parameter space for the non-Gaussianities in single-fleld inflation the shape is very 
different from the equilateral one. 

We flnd that the cosine with the local shape is also very small in the same region (63 ~ —5) 
where the shape is different from equilateral. As a result both the analysis of and /^^''' 
that have been carried on so far, have been largely insensitive to this region of parameter space 
and therefore there could be an undetected large signal in those triangular conflgurations. 
Notice also that in Fig. 1 the shape for which the scalar product in (23) is exactly orthogonal 
to the equilateral shape is peaked on both equilateral and on flat triangles, with opposite sign. 
This ensures that in the scalar product with the equilateral shape, there is a cancellation and 
the result is zero. Though the region where the shape of the non-Gaussianity is very different 
from the equilateral one is not large, being due to a partial cancellation of the two shapes 
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Figure 1: The shape of single-field inflation. Top Left: F^(^Q.^y (corresponding to C3 = 0), which 
is very similar to the template Equilateral shape. Top Right: Orthogonal shape: C3 = —5.4. The 
cosine of this shape with the equilateral shape is approximately zero. Bottom Left: Flat shape: 
C3 = —6. This shape is peaked on flat triangles where the two smallest fc's are equal to half the 
larger one, instead of on equilateral triangles. Bottom Right: F^s, which correponds to the case 
1 ^ Ical ^ 0(10): the contribution on flat triangles is clearly larger than in the case of F^(g.^)2. 



Cos( F,f,,,^,. ) Cos (F, 

1.0. 1.0 




Figure 2: Left: Cosine of single- field shape with the equilateral shape as we vary C3 with Cg <^ 1, the 
regime in which it is independent of c^. The two horizontal lines represent when the scalar product 
is equal to ±0.7, to give a rough measure of when the cosine becomes small. Right: Cosine with the 
local shape. 
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F local 

1.0 




Figure 3: The local shape. 



Fir{diwy and F^s, we consider it to be still roughly an fraction of the parameter space 
which deserves to be explored. 

One general characteristic of the shape of non-Gaussianities generated by a single-field 
inflation is that unless a large deviation from scale invariance is detected in the power spectrum 
the signal is always small in the squeezed triangle limit, where the lowest of the /c's is much 
smaller than the other two. This is the consequence of a theorem that says that the signal 
in that limit is proportional to the tilt of the power spectrum [1, 9, 10]. We stress that 
the theorem of [1, 9, 10] does not apply to inflationary models where fields other than the 
inflaton play a relevant role, either by producing entropy perturbations [6, 7, 16], or by 
slowing down the field with interactions [21, 28]. These models are allowed to produce large 
non-Gaussianties whose shape is peaked on squeezed configurations (like the one in Fig. 3), 
but they can also produce a signal peaked on equilateral or even flat configurations, as in [21]. 

In [26, 29], it was argued that modifications of the initial state of the inflaton field fluctu- 
ations can induce a departure from Gaussianity, and the shape of the induced signal, whose 
size strongly depends on the cutoff, was computed. Though we find the theoretical motiva- 
tion for these models quite unclear, the resulting shape of the signal is very similar to the one 
produced by single field inflation with a particular value of C3. Therefore, as we will explain 
in the next section, our analysis will automatically cover the signal expected in this case as 
well. 

Summarizing, the general analysis using the effective field theory of inflation shows us 
that at least when the system is away from de Sitter the non-Gaussian signal generated by 
single-field inflation when an approximate shift symmetry protects the Goldstone boson is a 
linear combination of two independent shapes. In general the resulting signal can be very 
different from the one associated with the equilateral shape. It therefore requires a more 
complete analysis that we will soon undertake. As we are now going to see, this is the case 
also when the system is close to de Sitter space. 
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2.1 Near-de- Sitter limit 



When the inequaUty (10) is satisfied, then the kinetic terms in the last fine of (8) become 
the dominant ones at horizon crossing. In this case, it is a good and useful approximation to 
take the exact de Sitter limit, sending H and Cg in (8) to zero in such a way that the quantity 
M|jij(l — Cg)/Cg approaches a constant value 

™Pi(l _ d) ^ 2M' . (25) 



C2 



Notice that the new parameter M is the same as defined in [10]. In this regime, the standard 
spatial kinetic term of the first line of (8) goes to zero, while the ones on the third line become 
important. The two kinetic operators possible in this regime are important in different regimes 
depending on the size of the ratio 4(^2 + d^Y^'^ /\di\. We will concentrate on the cases where 
only one of the two operators dominate at horizon crossing, neglecting the quite tuned case 
where both of the operators contribute at the same time, since, as it will become clear later, 
in this case we do not expect any new qualitative feature in the shape of the resulting non- 
Gaussianity. 

In the limit when \d\\ < 4(c?2 + d^Y^"^, then it is the higher-derivative spatial-kinetic term 
that dominates at horizon crossing. This is the regime of Ghost inflation [3, 10, 18, 19], that 
here we are able to describe as a continuous deformation from more standard models. In this 
case, the dispersion relation of the tt mode becomes extremely non-relativistic: 

Given the non-linear dispersion relation, the way an operator scales with energy does not 
coincide with its mass dimension as in the Lorentz invariant case. A rescaling of the energy 
by a factor s, E ^ s E, (equivalent to a time rescaling t — )> s~^t), must go together with 
a momentum transformation k — )• s^^'^k (x — )• s~^^'^x on the spatial coordinates) and a tt 
transformation tt s^^^w. Since in making predictions for inflation, we are interested in 
energy scales of order H, there are two interaction operators which are the most relevant 
(or, in technical language, the least irrelevant): 7r(9i7r)^ and (9|7r)(9i7r)^, which have scahng 
dimension equal to 1/4. For example, the operator ir^ we considered before has scaling 
dimension 5/4 and it therefore gives rise to an effect much smaller than the one induced 
by 7r(9i7r)^ and ((9|7r)((9j7r)^. The form of the non-Gaussianity these two operators induce 
is rather complicated, as the non-linear dispersion relation of (26) makes the wavefunction 
of TT, and therefore the resulting shape of the induced non-Guassianity, very complicated 
and not-writable in a closed form. The calculation can be easily carried on following the 
steps presented in [3, 18, 19], where the non-Gaussianity induced by the operator n{diny is 
computed. The wavefunction for tt is given by [3, 18]: 



-^ir) = -J^^i-rrH^>, I *^r= | . (27) 
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where r is the conformal time, and Hi^^ is the Hankel function of the first kind. Here 
we have expanded the operator tt in creation and annihilation operators in the usual way 
'^k — T^k^k + T^-k* ■ '^^^ contribution to the 3-point function from the operator ^{diny is 
given by the usual expression 



(28) 



■ / :^^''^i*(^)''S*(^X*(^)(^2 • ^3) + permutations + c.c. , 

where the sum above includes all symmetric permutations of the three momenta and the 
contour of integration should be rotated into the complex plane to ensure convergence as 
r — )■ —00. C^(a^7r)2 is the coefficient in the Lagrangian (8) of the operator 7r(5i7r)^/a^, which 
is equal in this case to 

^(9,.)^ = . (29) 

Analogously for the operator {df'!r){dj'!r)'^ we have: 

i'^A'^ks)(df^)(o,.r ^ (2^)''^^'^ {Y. ^(9f.Kd,.y (30) 

= -i (27r)^5(^) k}j 0) ' i/^%.)(a,.)^<(0)<(0)<(0) 

• / dr kl'Kl[*{T)7rl[*{T)7rl[*{T){k2 ■ h) + permutations + c.c. , 

J —00 

where C(^Q2^-^^g^^y is in this case given by 

The total 3-point function is given in this case by 

i 

Unfortunately, due to the complicated form of the wavef unction in (27), there is no closed 
form expression for the above integrals, which need to be integrated numerically. What is 
interesting to us for what concerns the data analysis, is that the resulting two independent 
shapes F^q^^-j2 and F(^q-2^-^(^q.^-^2 are peaked on equilateral configuration, but still different. 
This means that the two-parameter space of non-Gaussianities they generate leads to a signal 
that, in some region, is very different from the equilateral kind, peaked on flat-triangular 
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configurations. The values of the /jvl parameters they induce, defined in the usual way, are 
given by: 

where we have used the normalization of the power spectrum 

*~ 25r(i/4)2 ■ (^2 + 4)2/4 ■ 



(34) 



These expressions are valid when \di\ < 4(4 + 4)^^^ and the inequality in (10) is satisfied. 
Notice that they imply that the speed of sound of the fluctuations, meant as the coefficient of 
the linear term in k in the dispertion relation at energies higher than H, should be bounded 

by 

Id < 10-^(4 + 4)'/'. (35) 

where the absolute value is used since for di < 0, cj. is negative. 

When instead di > 4(4 + 4)^^^) then the spatial kinetic term at horizon crossing is a 
standard two derivative one, giving a speed of sound for the fluctuations equal to [3, 10, 20]: 

We will deal here with the case of positive di, and come back to the case of di < in the 
next subsection. Because of the extremely low speed of sound, the higher-derivative tri-linear 
operator of the last line of (8) becomes as important as the ones we considered in the former 
section. The shape that is generated by this operator is this time writable in closed form, 
and it has the following form 

5 1 

^(a?7r)(ai7r)2(^l, /cs) = • A| (37) 

(24 Ks^ - 4 K2^K^^K^ - 4 K^'^K^^ + 11 K^^K-^^^ - 3 Xs^i^i^ + K^^) 



Notice that, quite surprisingly, this shape is a linear combination of the shapes F^g.^^-^^ and 



F^s: 



-^(c)|7r)(<9,7r)2 



5 1 
3^ 



■Ali-^^F 



54 



7r(c)i7r)2 



81" 



fNL- Al=l 



(38) 



The non-gaussianity induced by the other two operators 7r(9.j7r)^ and tt^ is as in eq. (16), 
which leads us to conclude that in this case the 3-point function is given by 



(39) 
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with the values of the /atl's, defined in the usual way, given by: 

= ■ ^ - -2-662 ■ ■ ^ . (40) 

/;l = -^f7X^ + ^)- -4.115 -10- 



243 \(d2 + d3)^ 2 J ■ {d2 + ds) 



162 c2 - • df' ■ 

where we have used the normalization of the power spectrum, which in this limit has the 
peculiar form of 

4608 

^» = ^4 ■ 

A notable difference with respect to (16) is the coefficient of the operator tt^. This is due to 
the fact that at energies higher than the crossing energy E^r = diH/^d^ + (is)^^^, the spatial 
kinetic term is dominated by the A;^ term. This means that the scaling dimensions of the 
operators tt^ and 7r(c}j7r)^ become different, with the operator -k^ having scaling dimension 
equal to 5/4 while 'k{diTTY has scaling dimension equal to 1/4. Since the cutofi^ of the theory 
in this case is given by A ~ M[d2 + d^Y^"^^ it is quite straightforward to see that the operator 
TT^ becomes strongly coupled at the same scale as 7r((9j7r)^ if we substitute 

cl ci\K) {d2 + d,Y' ^ ^ 

with 63 of order one. Unless ^2 + d^ is small, the non-Gaussianity induced by this operator is 
rather negligible. 

We have three operators generating non-Gaussianities, but only two independent coeffi- 
cients (which means we have only a bi-dimensional space of non-Gaussianity). This can be 
clearly realized by writing the Lagrangian (8) in terms of c^: di disappears, and one is left 
with an expression of the form (neglecting numerical coefficients) 

M\i:' - clid^nf) + M^7r(9,7r^) + M^^^-|^7r3 + M^|(9j7r)(9,7r2) . (43) 

Since non-Gaussianities are generated when the modes cross the horizon, a time derivative 
contributes as a factor of while a spatial derivative as H/cg. Taking this into account, it 
is immediate to see that the operators 7r(9j7r)^ and {d'jTi){diTiY give rise to an f]\jL which is 
parametrically the same (^(l/c^). Notice again that the above expressions for /atl's are valid 
only when di > 4(^2 + c^s)^^^ and the inequality (10) is satisfied. 

Summarizing, the induced non-Gaussianity when we are close to the de Sitter limit can 
be large and detectable. The effective field theory shows that, depending on which operator 
dominates at horizon crossing there are two independent shapes of non-Gaussianities. As we 
concluded in the case where the system is less close to de Sitter, a more complete analysis is 
required. 



15 



2.2 Negative c 



There is one last regime which needs to be explored: the case where the squared speed of 
sound of the fluctuations is negative. The speed of sound is given by: 

where here we have not assumed that H is negligible. Imposing the time-kinetic term to be 
positive forces the inequality 2M^ > HM^^ + 3diHM^/2. We see that cj. can be negative 
when H is positive or di is negative. In this case, the dispersion relation, deep inside the 
horizon, is given by 

2 

PI 



where in the last term we have neglected a generally small correction proportional to HMj 
and to diHM^. When this happens, and the term in c^k'^ begins to dominate over the term 
proportional to k'^ at horizon crossing, the system is unstable, and the modes begin to grow 
exponentially. This does not mean the inflationary model is inconsistent. In fact, if the term 
in is large enough to dominate deep inside the horizon before the cutoff scale, the modes 
are stable in the ultraviolet, and then, as they redshift down to the Hubble scale, they become 
unstable and begin to grow exponentially. However, they do so only in a window of energies 
from when the term in c^k^ dominates down to the Hubble scale, when the modes freeze out 
and stop to grow. This means that such inflationary models are consistent. One interesting 
feature of these models is that they allow, if H is positive, to have a consistent violation of 
the null energy condition, and an inflationary model with a potentially detectable blue tilt of 
gravity waves [19, 20] 

Notice however that, as we will see, the fact that the modes grow exponentially for a 
window of time makes the induced non-Gaussianities rather large, and a consistent fraction 
of the parameter space of these models, as we will verify in the next section, is ruled out. 
This is why the operator in k'^ has to dominate soon enough at high energies, which explains 
why we still need to be close to the de Sitter regime, with the inequahty (10) not violated by 
too much. 

We will concentrate on two cases separately: when the \di\ > ~j^3h^j with di < 0, then 
the speed of sound is approximately given by 

c;.M<o, (46) 

and then in the case when |(ii| ^ ~^^^-> with H > 0, in which case 

e2^_™fi<0. (47) 



^Notice that H can be positive but small enough so that there is no exponential growth of the modes 
before horizon crossing [20] . This model would still imply a potentially measurable blue tilt of gravity waves. 
For what concerns the power spectrum and the non-Guassianities, this model is included in what studied in 
sec. 2.1. 
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Here we have assumed Cg <^ 1 and diH <^ M for simplicity. In the first case, as we saw 
in the former subsection, the condition for the term in c^/c^ to dominate over the one in 
imphes di < —4 (d2 + c?3)^/^, where, as usual, the 0{1) coefficients are not under control. The 
calculation for the power spectrum and the non-Gaussianities follows very closely the one of 
[19], and we do not give the details here. Deep inside the horizon, the dispersion relation 
is dominated by the stable term, which defines a stable vacuum, and the wave function 
agrees with the one of Ghost inflation (27). As far as the non-Gaussianities are concerned, 
since the modes at horizon crossing are dominated by a linear dispersion relation, the shape of 
the non-Gaussianities is exactly equal to the one with a positive given in eq. (22) and (37), 
with different values for the /atl's. These are given by: 

fif''^' = ^^^"^^'^^ 2.528 -10-^^^, (48) 
^^m^r = :^Ae"^^^ = 3.867. 10-- 



162|c,|2" |c.J7 



where in the second passage we have used that the normalization of the power spectrum is 
given by: 



4608 



|5 



A$ = 2{d2+d3)^/2 . (49) 

Notice, as expected, the exponential dependence on the ratio —di/{d2 + d^y^"^, with di < 0, 
that controls the amount of exponential grow of the modes before horizon crossing. Here we 
have assumed that the contribution from the operator tt^ is irrelevant. This is justified by 

the fact that, as we have seen in eq. (42), its expected importance decreases as we decrease 
the interval in energies between when the k^ term dominates and when the modes crosses 
the horizon. This is the interval during which the term in c^A;^ dominates. Because of the 
exponential dependence of /^r^ on this same interval of energies, we expect (and we will later 
verify) that this interval has to be rather small, and that therefore the operator tt^ gives a 
negligible contribution. 

The situation in the case where < "m^h^' "^i^^ > 0, is very similar. The condition 
for the c^/c^ term to dominate at horizon crossing becomes (^2 + ds) < SM^ / {H'^ H Mp^) . For 
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the non-Gaussianities, we obtain: 



J NL 



(a|7r)(9i7r)2 
JNL 



85 1 
324^ 



65 1 



Exp 



27/4 



1/4 



—di 



162 |cJ2 / 4ijMg, \ 
M^H ) 

-1.177-10-8 



Exp 



27/4 



15/2 



= -1.177-10" 



1 3/2 



1/2 



3/4 



13/2 



1/4 



5.179-10" 



1 3/2 



d 



1/2 



-1/4 



(50) 



where similarly in the second passages we have used that the normalization of the power 
spectrum is given by: 



9 



100 HM^,\c, 



Exp 



27/4 



1/4 



13/2 



1/2 



(51) 



and we have similarly neglected the contribution from the operator tt^. 

Summarizing, we see that there are consistent inflationary models with a negative squared 
speed of sound at horizon crossing. They induce exponentially large non-Gaussianities with, 
in general, two independent shapes. As for the other regimes that we discussed in the two 
former sections, a more complete analysis is required. 



3 Templates for single-field inflation non-Gaussianties 

In the former section we have studied the non-Gaussianities that can be generated by single 
field inflation, under the assumption that an approximate shift symmetry protects the Gold- 
stone boson. We have seen that the interesting regimes can be differentiated according to 
wether the system is close or not to dc Sitter, following the inequality in (10). We will show 
however that the same data analysis technique is sufficient for all these cases. 

Let us start from the case where the system is not close to de Sitter. In this case we have 
seen that the non-Gaussianities are given by a linear combination of the signal induced by 
two operators: tt^ and 7r(9i7r)2. Each of those two operators gives rise to shapes for the non- 
Gaussianities that arc peaked on equilateral configurations but are still quite different. Since 
non-Gaussianities from a generic single-field model of infiation which is not very close to de 
Sitter will be the result of a combination of the effect from these two operators, the resulting 
shape of the signal can be very different, peaked on fiat configurations or even changing sign 
as we go from equilateral to fiat triangular configurations. We therefore conclude that it is 
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necessary to do the analysis of the non-Gaussianity jointly for both of the shapes Fj^(^Q.^Y and 
F^3 . The numerical analysis becomes much simpler if instead of doing the analysis for precisely 
these two shapes, we use templates which are very similar but which, being factorizable as 
a product of functions of fci, k2 and ^3 or as a sum of a small number of terms with this 
property, are computationally efficient [30]. For the equilateral shape, this is given by [12]: 

-^equii.(A:i, A;2, h) = f^^'^' ■ 6A| • ( -pp - pp - pp - ppp + k p7,3 + perm. 

\ fv-|^n'2 rt-j^/i'^ fV2'i'3 nj'^i\"2'^^ '^1"'2 3 

(52) 

where the permutations act only on the last term in parenthesis. The second independent 
template we choose is similar to the orthogonal shape (see upper-right panel of Fig. (1)), and 
it is given by: 

Fo.thog.{h,k2,k3) = /°'1'^°^ -6A|- ( -T^s - pis - - ppp + ppp + (5 perm.) J , 

(53) 

where the permutations act only on the term immediately to the left. In practice, Fequii. is 
a good template for the equilateral shape with 63 ~ 0, while Forthog. is a good template for 
the orthogonal shape, 63 ~ —5.4. As we will later see, these two templates will be able to 
reproduce with good accuracy the bidimensional space spanned by the exact single-field shape 
as we vary the contribution of F^(g.„)2 and by varying and 63. 

Notice that the orthogonality of the shapes is in reality experiment dependent: we can 
consider not only the three dimensional cosine but also in the case of a CMB experiment, the 
2D cosine between the bispectra Bi^i^i^ as [2]: 

where we have defined the 2D scalar product: 
with fi^ 

,h,h being a combinatorial factor equal to 1 if the three Vs are different, to 2 if two of 
them are equal and to 6 if all of them are equal, and with Zmax being the maximum I of the 
CMB survey. Here the bispectrum is defined in such a way as 

{O'hmiO'hm^f^hrm) = ( ^ ^ ^ ) Bi^i^i^ , (56) 
\ nil ^2 ^3 / 

while 

{(^Iimia*l2m2) = C'h^/ii2^mim2 , (57) 

where the a^^'s are the CMB multipoles. The 2D cosine is expected to give a quantitative 
measure of the correlation of two shapes after the transfer functions and the sky projection 
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are applied. If we define a more general family of orthogonal shapes 

1 



orthog 



orthog. 
NL 



1 



-Foquil.(fcl, h 



-6 



(58) 



we can consider not only the 3D cosine of eq. (24), but also the /max-dependent 2D cosine, 
as we let the constant c in eq. (58) vary. The value of the various cosines between Forthog. 
and -Fequii. as we vary c is shown in Fig. 4, where one can see that the value changes as we 
pass from 3D to 2D, and as we let /max vary. We choose c = 2/3 in eq. (58) because this 
represents a good approximation to a 2D orthogonal shape for an experiment like Planck. 
We stress that, if one takes into account of the correlation between Forthog. and -Fgquii., as we 
will do, all of the different choices of c are equivalent, because they represent different linear 
combinations of the same shapes. 
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Figure 4: Cosine between Fgquii. and Forthog. as we let c vary and as we use the 3D and 2D cosine, 
or as we vary the /max in the definition of the 2D cosine. As expectable, we see that how much the 
two shapes are similar depends on the survey. Our choice of the template Forthog. corresponds to 
c = 2/3, which is quite close to a 2D orthogonal template with repect to Fequii. at Planck resolution. 

We find that the template built from Fgquii. and -Forthog. for estimating the full one- 
parameter family of shapes in single-field inflation is: 

Ftemplate(/i^l; k^, Cg, C3) = /^^^^'(c^, C3) F^qnilXklj ^2; k^) + f'NL^^^i'^s, C3)-Forthog. (/i^l ; ^^^2? k^) . 

(59) 

Here /^^'''(cs, C3) and /^^^^^"^ (^s, C3) are chosen by assuming that the expectation values of 
the estimators for our factorizable templates are given by the Fisher matrix prediction using 
the 3D scalar product defined in (23): 



20 



' \ \ ^orthog. '^orthog. / \ -^orthog. orthog. / / f ,^2 ^ 



1.040 1.210 
0.1079 -0.06572 




(60) 



where f^^^^^ (csjCa) and f^Li^s^Cz) are given by eq. (21). This choice is dictated by the 
definition for the estimators for /^"'^' and f^^^^^' that we will make in the next section. See 
App. A for a more straightforward but still completely equivalent alternative definition. 

Quite remarkably, the cosine of this template with the exact shape for any value of and 
C3 is always larger than 0.91. Notice that this cosine is defined using the /^^''' and 
from eq. (117) of App. A. This is in fact the definition of /^^''' and f'^^i°^' that minimizes the 
norm of the shape F — ^template with respect to the 3D scalar product. This is in fact the choice 
that makes the template shape the best approximation to the one of single field inflation for 
any value of Cs and 63 ^. The cosine is shown in Fig. 5, where we see that it is always very close 
to one. This tells us that our two independent shapes are covering the parameter space of the 
shapes generated in single-field infiation quite well. The scalar product with the exact single- 
field shape is minimal when 63 ~ —5.4, where it is approximately equal to 0.91. This value 
of C3 corresponds to when the shape of the exact single-field non-Gaussianities approaches 
the orthogonal shape (see upper-right panel of Fig. 1). This means that our template of the 
orthogonal shape is not extremely accurate, though it is still a very good approximation. 
Errors up to 9% are satisfactory in the absence of a detection of a non-Gaussian signal. The 
main reason why our approximation is not better is that the template Forthog.(fci, /^2, /^s) does 
not have the correct behavior in the squeezed limit where ki <^ A;2, /C3. In this regime, the 
single-field shape F{k-i, k2,ks) is of order 0{l/ki), while our orthogonal template is of order 
0{l/kf) (see Fig. 6 and 7, where again we use the /^"''' and f^^j^*^^' from eq. (117)). This is 
not a major problem, as what is the relevant number for the similarity of the shapes is the 
value of the cosine. In App. B we show a simple generalization of the template we use, which 
has a cosine with the exact shape always larger that 0.99. Though clearly extremely good, 
this generalized template is rather difficult to implement numerically and therefore given that 
our template is already a good approximation to F in all the parameter space and that we 
will have no detection of non-Gaussianity we believe our analysis is sufficient. 

By inverting the above relationship (60), we obtain the estimate for and £3 in terms of 
fNL^ and f^L^"^ - These expressions will allow us to translate the limits on the /nl parameters 
into limits directly on the parameters of the Goldstone boson Lagrangian (8) ^. 

^While this choice of /^^''' and f^^"^' is better for this kind of comparison, it is just a hnear transforma- 
tion of the /tv^'s defined in eq. (60) that affects neither the minimum of the cosine, nor, as shown in App. A, 
the constraints on the parameters of the Lagrangian. 

^In [31], it was argued that the shape of the non-Gaussianity induced by the initial state of each mode 
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Figure 5: The cosine between the template shape ^template and the exact single- field shape F as we 
vary C3, for Cg ^ 1 where it is independent of c^. The cosine is always very close to one, reaching its 
minimum, equal to approximately 0.91, for C3 ~ —5.5. To help visualization, we plot also the line 
corresponding to a cosine equal to 0.9. 
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Figure 6: Left: ^template for £3 = 0. This plot should be compared with the left-top panel of Fig. (1). 
Right: ^template for C3 = — 5. This plot should be compared with the right-top panel of Fig. (1). 
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Figure 7: Difference between the template shape ^template and the exact shape F for £3 = (on the 
left) and for C3 = —5 (on the right). Notice that the two figures have different scales on the vertical 
axes. 
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3.1 Near-de- Sitter limit 



When the inflationary spacctimc happens to be very close to de Sitter, the inequality in (10) 
can be satisfied, and the behavior of the system changes. In sec. 2.1, we have seen than 
there are two possible regimes: one in which the goldstone boson has a linear dispersion 
relation uj — Cgk, with a small speed of sound [3, 10, 20], and one in which the dispersion 
relation is quadratic in the wavenumber uj oa k^. This last regime is known as Ghost Inflation 
[3, 10, 18, 19]. If di > 4:{d2+d3y^^, then the operator which induces a linear dispersion relation 
dominates. In this case, we have seen that there are three operators which can induce large 
non-Gaussianities. The resulting non-Gaussianity is dependent on just two combinations 
of coefficients: the speed of sound Cg (equivalent to the parameter di) and the coefficient 
Cz/{d2 + (is)'^. The same discussion of the former section applies also in this case: the shape 
induced by each of the operators is peaked on equilateral configurations, but the total signal 
can in principle be peaked on even flat-triangle configuration. This requires a more complete 
analysis that is able to constrain the signal peaked on fiat-triangles. To this purpose, the two 
templates we defined in the former section are enough and we just need to give the Fisher 
matrix for this case as well: 



/■orthog. 

Jnl 



idl,C3) 



Y -^equil. '-^equil. J y -^equil. '-^equil. y 



equil. 



^■7T(a-7r)2 "-^orthog. 
^orthog. '-^orthog. 



irS "-^orthog 



- orthog. '-'^orthog. 



1.040 1.210 0.9878 
0.1079 -0.06572 0.1613 



equil. '-^ equil. 

-^(a27r)(07r)2 "-^orthog. 
^orthog. "^orthog. 



fNLih/id^ + d,Y) 

(a2^)(9,7r)2 



Jnl 

f%\{c3/{d2 + d3 



(62) 



\ f. 



NL 



(di) J 



\ f. 



NL 



(di) j 



This relationship is valid when the inequality (10) is satisfied and d\ > 4(^2 + d^]^!"^ . Notice 
that since there are two independent parameters, d\ and £3/(^2 + (^3)^, the non-Guassianity 
in this case can be written as a linear combination of two shapes, each one proportional to 
l/(i^^^ and £3/(^2 + c^s)^- The resulting Fisher matrix can be simply obtained from the one 

as not being in the standard Bunch-Davies vacuum is well approximated by what the authors of [31] call 
the Enfolded template. Though we find the theoretical motivation of this effect quite unclear, our parameter 
space is so general that it includes the enfolded shape, which is given by the following choice of the /jvl 
parameters: 

fequil. 1 j^enf. j^orthog. 1 j-enf. ( a-\ \ 

Jnl — Jnl > Jnl ~ ~i2 ■'^^ ' ^ ^ 

where we have defined /™| A|, to be the coefficient that multiplies the enfolded shape presented in [31]. 

^°Notice that for the particular sub-case we arc talking about, since the shape Ff^Q.^^^^^Q.^^^i is a linear 
combination of the other two shapes Fj^(^Q.T^Y and F^3, we can use the result of the former section and 
conclude that the cosine of our template with the one produced by these models of inflation will be larger 
than 0.91. 
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above after using eq. (40) , and it reads 



/■orthog. 
JNL 



((^1,63) 



-6.791 • 10=^ 
-9.441 • 10^ 



-4.979 • 10-2 
2.705 • 10-3 



(63) 



Still close to dc Sitter, but in the opposite limit where \di\ < 4(^2 + ^3)^^^, the system 
approaches Ghost inflation. In this case the same discussion of the case di > 4(^2 + d^Y^'^ 
applies, with the only difference that now the operator tt^ is always irrelevant. In particular 
the two templates we use are still a good approximation: the cosine never goes below 0.9. 
The resulting Fisher matrix is given by 
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(64) 



Notice that the bar over F, 



7r(97r)2 



and F(g2^)(g^)2 



means that those are the shapes generated 
in Ghost inflation respectively by the operators ^{dny and (9|7r)(97r)2. Because of the 
peculiar dispersion relation ou oa k^, the form of this shapes is different from the ones of 
eq. (22) and (37), and can only be computed numerically as in eq. (28) and (30). We 
remind the reader that the above relationship is valid when the inequality (10) is satisfled 
and \d,\<A{d2 + d3y/\ 



3.2 Negative cj 

In sec. 2.2 we saw that still quite close to the de Sitter hmit, when if > or when di < 0, 
the squared speed of sound of the fluctuations at horizon crossing can be negative. In this 
case the non-Gaussianities tend to be exponentially large. In these models the inequality (10) 
can be violated, but not too much, as the operator in has to dominate quite early at high 
energies in order not to have an exponentially large level of non-Gaussianities. In this sense, 
also these models are quite close to de Sitter space. 

In the case in which di is negative, with di < —4(^2 + d^)^ and \di\ > 4\H\Mp^/{M^H), 
then the operator that dominates at horizon crossing is the one proportional to di. In this 
case, there are two shapes for the non-Gaussianities that can be large, and the Fisher matrix 
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relation with our templates is given by: 



/•orthog. 
JNL 



{dl, \Cs 



^■k(dTT)'^ '-Pcquil. 



^cquil. "-^cquil 

/ -^7r(97r)2 --^orthog. \ 
\ -^orthog. '-^orthog. / 

1.040 0.9878 
0.1079 0.1613 



-^(a^7r)(S7r)2'-^cquil. 
^cquil. "-^equil. 
(a|7r)(a7r)2 
-^orthog. '-^ortho 



vr(9i7r)2 
NL 



NL 



J 



{dl, \cs\) 
{dl, \cs\) 



lft^\d„\cs\) 

NL 



{dl, \cs) 



(65) 



where fpl^'""^ 



(rfi, |cs|) and fj^l ' (c^IjIcsI) are given by eq. (48). It turns out that the 
correlation of the non-Gaussianity generated in this particular model with the equilateral 
template is always very large (the cosine is greater than 0.999). For this reason, we will 
perform the analysis for this model using only the equilateral template. 

The other case in which it is possible to have < is when H > 0, \di\ < AHM^J (M^H), 
and ((^2 + ^3) < SM^ / {H'^ H Mp^) . In this case, the kinetic operator that dominates at horizon 
crossing is the one proportional to H, and the Fisher matrix is equal to the one above, with 
the replacement of the /atl's with the one of eq. (50). Notice that this inflationary model 
predicts a blue tilt of gravity waves [19, 20]. 

Summarizing, Eqs. (60), (62), (64) and (65) are what is necessary to translate the limits we 
will obtain on the templates Fcquu. and Forthog. into limits on the parameters of the Lagrangian 
of single field infiation (8), as we are now going to do. 



4 Results from WMAP 
4.1 Analysis pipeline 

We use the pipeline from [11] to obtain constraints on /^^'^' and f^^j^"^' from the 5-year 
WMAP data. Since the relevant bispectra are factorizable in ki, k2, k^, the optimal estimators 
can be written down in a straightforward way following [30, 12, 32]. We define the functions: 

= J'^-^Aj{k)je{kr) (66) 

-J^Aj{k)P{k)je{kr) (67) 

le{r) = P-^Ajik)Pik)y%ikr) (68) 

Se{r) = P-^Aj{k)P{kme{kr) (69) 
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A{n,r) 


= ^a,(r)(C- 






em 




B{n,r) 




^a)e,rnyem{^) 




em 




C(n, r) 












L>(n,r) 







where P{k) is the initial power spectrum of the Newtonian potential $ and Aj(A;) is the 
temperature transfer function. 

We analyze the six five-year WMAP maps in V+W channels at Healpix resolution A^side = 
512, using the kq75 mask described in [33]. These maps are reduced to a single harmonic- 
space map {C~^a)irn as described in Appendix A of [11]. The inverse (signal + noise) filter 
= {S + N)~^ includes the sky cut, inhomogenous noise, and channel-dependent beam 
in the definition of the noise covariance. We also marginalize the monopolc, dipole and 
foreground templates for synchrotron, free-free and dust independently in each of the six 
channels. 

For each value of r, we define pixel-space maps A,B,C,D by filtering the WMAP data 
(C~^a) as follows: 

(70) 
(71) 
(72) 

(73) 

em 

The optimal estimators f^^^' and f^^°^' are then given by: 

Inl'- = ^ /drr^[-3A(n,r)S(n,r)2 + 6S(n,r)(7(n,r)L>(n,r)-2i^(n,r)^ 

+6 {A{n, r)B{n, r))^^ B{n, r) + 3 {B{n, r)B{n, r))^^ A{n, r) 
-6 (S(n, r)C7(n, r))^,^ ^(n, r) - 6 {B{n, r)D{n, r))^^ C{n, r) 
-6 {B{n, r)C{n, r))^^ D{n, r) + 6 {D{n, r)D{n, r))^^ D{n, r)] (74) 

fNL°^- = — ^ /cirr2[-9A(n,r)S(n,r)2 + 18S(n,r)C'(n,r)L'(n,r)-8L>(n,r)^ 

-''orthog J 

+18 {A{n, r)B{n, r))^,^ B{n, r) + 9 {B{n, r)B{% r))^,^ A(n, r) 

-18 (fi(fi, r)C(fi, r)) D{% r) 18 (fi(fi, r)D(S, r))^,^ C(fi, r) 

-18 (i?(n, r)C(n, r))^,^ ^(n, r) + 24 (D(n, r)D(n, r)),^^ ^(n, (75) 

Each estimator contains a 3-point term and a 1-point term which improves the variance of the 
estimator in the presence of inhomogeneous noise or a sky cut, as described in [12]. Whenever 
an expression such as (^(n, r)S(n, r))MC appears in Eqs. (74), (75), it denotes a Monte Carlo 

average in which we construct the filtered maps (A(n, r) and S(n,r) in this case) from signal 
-|- noise simulations instead of the WMAP data, and take the average over many simulations. 

In implementation, the r integrals in the definitions of f'^^^^', f'^Ji^"^' must be replaced by 
finite sums. We use the optimization algorithm from [13] to minimize the number of sampling 
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points, finding that 30 and 54 r- values are necessary, for the equilateral and orthogonal shapes 
respectively, to preserve the Fisher matrix to one part in 10^. 

We compute the normalization constants A'equii, -/Vorthog by running the pipeline on ensem- 
bles of non-Gaussian simulations, using the simulation algorithm from [13]. This ensures that 
each estimator is an unbiased estimator of the corresponding /jvl parameter, assuming no 
additional contributions to the bispectrum 

{InT-) = InT- (assuming = = = 0) 
(Ar^-) = Cl"'- (assuming = /^7^- = ... = 0) (76) 

(We have written "• • • " here to indicate other possible terms in the 3-point function, for 
example point sources.) The covariance of the estimators f^l^^', f^T^"^ ^^^^ determined by 
Monte Carlo. 

There is an alternate version of this construction, in which the estimators f^i'^', fNL^°^ i 
■ ■ ■ are defined differently (by taking linear combinations) in such a way that each estimator 
has unit response to the corresponding /jvl parameter and zero response to the other (A^— 1) 
parameters. Details of this construction are given in Appendix A. For now we note in advance 
that the analysis of single-field infiation in the following subsection will not depend on which 
definition is used; it is simply an arbitrary choice that does not affect parameter constraints. 
We have chosen to use the definition in Eqs. (74), (75) for consistency with previous analyses 
[12, 47, 22] so that our hmits on f^l'^' in the next subsection can be directly compared to 
results in the literature. 



4.2 WMAP constraints on f^f', ff^""^- 

We encountered one issue in our analysis which affects previously reported constraints on 
/^^''■. In principle, the r integral in Eqs. (74), (75) should run from r = to r = oo. Previous 
analyses [12, 47] have cut off the integral at the horizon rhoriz ~ 14500 Mpc. Empirically, 
we find that setting rmax ~ '"horiz + (300 Mpc) is needed for convergence; truncating at the 
horizon significantly underestimates the error cr(/^^'^ ). 

The need to extend the r-integral beyond the horizon does not indicate that causality is 
violated. The r integral in the estimator arises from writing a delta function {2nY'6'^{ki + k2 + 
k^) as an integral J d^re^^'^. so formal contributions to the estimator from r > rhoriz do not 
correspond to physical contributions from outside the causal horizon. However, one can use 
causality to show that Q;^(r) = for r > rhoriz- This implies that the estimator for the local 
shape can safely be truncated at the horizon, so this issue does not affect previously reported 
estimates of f]^l (e.g. [47, 48, 22]). 

^^Some secondary contributions to Jnl have been studied and have all been predicted to be small compared 
to the WMAP statistical error [11, 13, 34, 35, 36, 37, 38, 39, 40, 41, 43, 44, 45, 46], although it is still not 
clear that all possible important secondaries have been studied. 
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optimal estimator 


suboptimal estimator 


WMAPl (kpO mask) 
WMAP3 (kpO mask) 
WMAP5 (kq75 mask) 


fZ'' = 125 ± 177 
/^7'- = 146 ± 149 
/^7^=155±140 


f^f- = 236 ± 204 
r^f- = 145 ± 162 



Table 1: Constraints on /^^''' from 1-year, 3-year and 5-year WMAP data. Previously reported 
results have underestimated the error o"(/^^''') by truncating the integral in the estimator (eq. (74)) 
at Tmax = '"horiz ^ 14500 Mpc, rather than taking rmax = c«. The results shown here for the optimal 
estimator represent the best possible statistical error for each dataset. We note that the statistical 
errors in WMAP3 and WMAP5 are roughly equal, even though WMAP5 is more sensitive, because 
the kq75 mask introduced in the 5-year analysis is more conservative than the kpO mask used in 
the 1-year and 3-year analyses. For comparison with previous results, we also show the constraint 
obtained using a suboptimal estimator similar to the one described in Appendix A of [22] . 

For the five-year WMAP data in (V+W) bands with the optimal estimator, kq75 mask 
and foreground template marginalization, we find /^^''' = 155 ±140 and f^]]^"^' = —149 ±110 
(errors are la), indicating no detection of non-Gaussianity using either shape. The 2-by-2 
covariance matrix is given by 



C = 



Var(/X') Cov(/^7\ 
Cov(/^t\/°r^-) Var(/°r^-) 

1.96 X 10^ 5.0 X 10^ 
5.0 X 10^ 1.21 X 10^ 



WMAP 

(77) 



where the subscript wmap reminds us that these quantities depend on the WMAP data. We 
find that the correlation Corr(/^^''', f'^^°^') in the Monte Carlo simulations is consistent with 
the ~ 0.32 cross correlation expected from computing the 2D cosine defined in eq. (54). 

Since our estimators are optimal, the la error o"(/^^'^') = 140 is the best constraint that 
can be obtained from the five-year WMAP data, if the r-integral in eq. (74) is properly taken 
to Tmax = oo. As a check, we verified that the statistical errors agree with the Fisher matrix 
(assuming f~^^'^ scaling). For comparison with previous results, we also report constraints on 
ff^f- from 1-year and 3-year WMAP data in Table 1. 



4.3 WMAP constraints: non near-de-Sitter models 

In the next few subsections, we will use the WMAP measurements f^f' = 155 ± 140, 
fNL^°^ ~ —149 ± 110 to place constraints on single-field inflation in different model spaces. 
In this subsection we will consider the parameter space (c^, 63) of the models not-so-near- 
de-Sitter described in sec. 2. In the following subsections we will analyze the near de Sitter 
models and the models with < 0. 

Our basic tool will be a statistic which quantifles agreement between the observed 
bispectrum and the model bispectrum. Given model parameters (0^,63), our estimators 
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/nl^^ fNL^°^ acquire expectation values given by: 




( 



10 
243 



1(1 

(l-l/cD(53 + 3c2 



/2) 



) 



(78) 



by combining Eqs. (21), (60) above. 

Given estimates (/^^^'Owmap, (/Ari,''°^')wMAP from the WMAP data, with the associated 
2-by-2 covariance matrix C of eq. (77), we define a statistic by: 



X (Cs, C3)wMAP 



v{Cs, C3)wMAP 



c 



1 



v{Cs, C3)wMAP 



(79) 



where: 



'^(Cii, C3)wMAP 




(/^r^-(c.,c3)) 




(80) 



We have used the notation x^(cs, C3)wmap to emphasize that the statistic depends both 
on the model (c^, 63) and the WMAP data, and measures agreement between the two. If 
we construct an ensemble of non-Gaussian simulations in a fixed model {cg, C3), and evaluate 
X^(cs, C3)sim on each simulation using the true (c^, C3) of the model, then it will be distributed 
as a random variable with two degrees of freedom (as suggested by the notation). Turning 
this around, we can test the hypothesis that a given (c^, £3) is consistent with the WMAP data, 
by comparing the value of x^(cs, C3)wmap to a distribution with two degrees of freedom. 

In Fig. 8, top panel, we show 68%, 95% and 99.7% confidence regions in the {cg, C3)-plane 
using the WMAP data. More precisely, we define the confidence region corresponding to 
p-value p (where p=0.32, 0.05, or 0.003) by the inequality 



where F is the CDF for a random variable with two degrees of freedom. It is seen that 
the simplest, non-interacting single-field model {cg = 1, 63 = 0) is consistent with the WMAP 
data (the p- value is p — 0.1). The allowed region of parameter space shows a nontrivial 
structure which arises because the two primordial bispectrum shapes F^(a^)2,F^3 are nearly 
degenerate. There is an allowed region in the lower left corner with C3 ^ —5.4 and small Cg. In 
this region, there is a near-cancellation between the bispectra associated with the operators 
7r(5i7r)^ and tt^ (see Fig. 2). In the limit Cg — )■ 1, the parameter £3 is not bounded, although 
the WMAP data do give the constraint -9920 < (1 - cDcs < 1000 (95% CL) in this hmit. 

^^Note that we define confidence regions in the sense of frequentist statistics: we say that a model (c^jCs) 
Ues outside the 95% confidence region if, in an ensemble of simulations with model parameters (0^,53), the 
value of x^(cs,C3)wMAP is larger than x^(cs)C3)sim 95% of the time. The p- value that we assign to a model 
(csjCs) should be interpreted as the significance of rejecting the null hypothesis that the model is consistent 
with the WMAP data, using as the test statistic, not as a probability of (05,03) being the true model. 
(For example, if we integrate p over the (cs,C3) plane, the integral will not be equal to one.) We have used 
hypothesis testing rather than Bayesian inference for this problem because there is a clear choice of test 
statistic x^j whereas there would be no clear choice of prior density P{cs, cs)dcsdcz in a Bayesian analysis. 



X^(Cs,C3)wMAP < F ^(1 -p) 



(81) 
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Figure 8: Top panel: 68%, 95% and 99.7% confidence regions in the single-field inflation param- 
eters (csjCs) from five-year WMAP data, obtained from an analysis which uses /^^'^' and 7^2^^°^ 
(eq. (81)). Bottom panel: Confidence regions obtained from an analysis using /^^'^' alone (eq. (82)), 
showing weaker constraints. 
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For comparison, in the bottom panel of Fig. 8, we show 68%, 95% and 99.7% confidence 
regions using an analysis which discards f^^^"^' and obtains constraints from f^l'^' alone. 
More precisely, we define an "equilateral" by: 



Var(/Z"-) 



,2/ ~\ K\JNL {"-3,^3)/ - KJnL ;WMAP; 

X {Cs: Cgjequil = . fcquil.x (^2) 



and define confidence regions by converting to a p-value using statistics with one degree 
of freedom. Using /^^'^' alone, it is seen that there is a degeneracy which allows models with 
63 —5.4 and arbitrarily small Cg. As we move along the degeneracy line in the — )■ 
direction, we get large contributions to /^^''' (proportional to (1 — l/c^) and c^/cj.) which 
nearly cancel. If f'^^j^'^^' is included in the analysis (top panel), the degeneracy is broken for 
sufficiently small Cg, because a detectably large f^^t"^ generated. 

We can also ask: what is the WMAP constraint on the sound speed Cg during infiation 
The answer is different depending on what assumptions we make about the parameter C3. If 
we assume 63 = 0, or DBI infiation (63 = 3(1 — c^)/2), we can discard f^^°^' and constrain Cg 
using /^^''' alone. This is because the bispectrum is always highly correlated to the equilateral 
shape, so that including f^^^"^' would not increase the statistical power. We find the following 
lower bounds: 

Cg > 0.048 (95% CL, C3 = 0) (83) 
Cg > 0.054 (95% CL, DBI infiation) (84) 

To get the most conservative constraint on Cg, we would marginalize C3 instead of assuming 
a specific form. More precisely, we define a one-variable statistic by evaluating the two- 
variable x^(cs) C3) at the value of C3 which minimizes the fo^' ^ given c^: 



X^(cs)marg = Hiiu (x^{cg, C3)wMAp) (£3 marginalized) (85) 

Converting this to a lower limit using statistics with one degree of freedom, we find: 

Cg > 0.011 (C3 marginalized, 95% CL) (86) 



This very general lower limit applies to all single-field infiation models regardless of the value 
of the coupling C3. To obtain it, it is necessary to include the new shape °^ ^^e analysis. 
Using f^^^^' alone, there is no lower limit on Cg when 63 is marginalized, due to the degeneracy 
seen in the bottom panel of Fig. 8. 

So far in this analysis we have neglected the metric fiuctuations. As shown in detail in 
[10], in these kind of models these give only subleading corrections suppressed by the slow 
roll parameters, and can be safely neglected. 

"^^We will constrain now for the models not-near-de-Sitter, and will come back again to Cg in the next 
subsection. 
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Figure 9: WMAP constraints in the (di, 03(^2 + ds) ^)-plane, for the near de Sitter model with 



4.4 WMAP constraints: near de Sitter models with di > 4:{d2+d^ 



1 1/2 



We can also use the two-template WMAP analysis to put constraints on the near de Sitter 
parameter spaces described in sec. 2.1. In this subsection, we will consider the case di > 

m + d,y/\ 

We take the parameters of the model to be ((ii, 63(^2 + ds)^'^). As in the preceding 
subsection, we define a statistic which quantifies agreement between the model parameters 
and the WMAP data, using eq. (63) above for the expectation values of the /nl estimators. 
In detail, is defined as a function of {di, £3(^2 + c^s)"^) by: 



(87) 



(Al"') ] - f -6-791 X 103 -4.979 x IQ-^ \ / df^^ 
ifNt'-) ) ~ V -9.441 xlO^ 2.705 x 10-3 J 53(^2 + 4) 

/ ?cquil.\ / ?equil.\ 

\Jnl I ~ Unl Jwmap 

; Porthog.\ / Porthog.\ 

\Jnl ) - Unl Jwmap 
= v^C-'v (where C = Cov(/^7'-,/°f°^-)) (89) 

In Fig. 9, we show la, 2a, and 3a contours in the {di,Cs{d2 + d^,)^^) plane obtained using this 
statistic. The non-interacting model {di ^ 1, C3{d2 + 6^3)"'^ = 0) is consistent with the 
data at 2a. (This is the case for all parameter spaces considered in this paper, since the of 
the non- interacting model is independent of the parameter space into which it is embedded.) 
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Figure 10: WMAP constraints in the {di, d2 + d^) plane, for the near de Sitter model in the ghost 
inflation limit \di\ < 4((i2 + ^3)^/^. Note that the range on the y-axis has been chosen so that only 
models which satisfy this inequality are shown. As described in the text, a Icr region does not appear 
in the plot because this parameter space does not include any model whose gives better than la 
agreement with the WMAP data. 

It is worth noticing that if we imposed d2 + (is to saturate the bound di > 4(^2 + ^3)^''^ and 
we also imposed 63 not to be too large, then many of the most extreme values of the y-axis 
in Fig. 9 would be excluded. Using eq. (41), the sound speed is given 
and lower limits on Cg are given by: 

Cs > 0.073 (95% CL, £3(^2 + 4)"^ = 0) (90) 
Cs > 0.0189 (95% CL, ^3(^2 + ^3)"^ marginalized) (91) 

Note that the requirement that < 1 implies di < 319 in this model. 

So far we have neglected metric fluctuations. Following [20, 10] it is straightforward to see 
that this is a good approximation for di < (Mpi/(30M))^. For the parameter space we are 
interested in, this translates in M < Mpi/10 to a good approximation. When this inequality 
is violated, as shown in [20, 10] the squared speed of sound of the fluctuations becomes 
negative and the modes begin to grow exponentially before crossing the horizon. In this case 
these models become similar to the ones shown in the next sec. 4.6, and, as we will see, are 
generically very disfavored by the data. 

4.5 Near-de-Sitter models with \di\ < 4(^2 + 0^3)^/^ 

The next space of models we consider is the near de Sitter limit with \di\ < 4(^2 + d^y^"^ 
relevant for Ghost inflation. In this case, we take the model parameters to be {di,d2 + ^3). 
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Combining Eqs. (33), (64) above, the expectation values of the f^L estimators are given by: 

{r^t\dud2 + ds)) \ / 0.8625 0.9685 \ / 138.1(^2 + a^s)"'/' \ ...^ 
(/^?°^ (di,d2 + d3)) J V 0.2621 0.1667 J [ nm{d,){d2 + d^)-''"' ) ^ ^ 

Using this expression, we can define a statistic as in Eqs. (89) above. In Fig. 10, we show 
2(7 and 3(T contours in the (^1,(^2 + dz) plane obtained using this statistic x^(di,d2 + dz). 
Note that no la contour appears in the figure. This is because there are no models in this 
parameter space which satisfy x^(di,d2 + ^^3) < 2.28. (This value of x^ corresponds to p- 
value p = 0.32, and is the threshhold value for the la region, as described in sec. 4.3.) For 
example, if we try to solve for the values of di, (^2 + d^) which satisfy x^ = (or equivalently, 
(/atT'') = (/atT^Owmap and {/'^l'"^') = (/^l''°^')wmap), then we find a negative value of 
{d2+d3)^^^ , so no model with x^ = exists. In contrast, the other parameter spaces considered 
in this paper do include models with x^ = 0, and a la region appears in the associated figures 
(Figs. 8, 9). The absence of a la region in Fig. 10 should not be interpreted as evidence that 
the near de Sitter model with < 4(^2 + dsY^"^ is inconsistent with or disfavored by the 
WMAP data, since it does include models which are consistent with WMAP at 2a. 

As in the former models, so far we have neglected metric fluctuations. Following again [20, 
10] it is straightforward to see that this is a good approximation for (0^2+^3) < {Mpi/ (30M))^°. 
For the parameter space we are interested in, this translates in M < Mpi/10 to a good 
approximation. Similarly to the former case, when this inequality is violated, as shown in 
[20, 10] the squared speed of sound of the fluctuations becomes negative and the modes 
begin to grow exponentially before crossing the horizon. In this case these models become 
again similar to the ones shown in the next sec. 4.6, and, as we will see, are generically very 
disfavored by the data. 

4.6 Near-de-Sitter models with < and -di > \4:HM^J{M^H)\ 

In the near de Sitter model with < and —di > {M^H)\, the bispectrum (eq. (48)) 

is always highly correlated to the equilateral shape. This is because the signs of f^^^^ and 
Inl^^^^^^ are always the same, so it is not possible to get the near-cancellation which leads to 
the orthogonal shape. (More precisely, we find that the cosine between the equilateral shape 
and the bispectrum in this model is equal to 0.9995 for small jcgj and becomes even closer to 
one as we let \cs\ grow.) For this reason, we will discard ffj^"^- and only use fj^f' in the 
analysis. Using eq. (65), and neglecting the orthogonal template, the expectation value of 
Inl^' is given by: 

{InT-) = ( 1-040 0.9878 ) J ,..3 (93) 

where here, contrary to what we did in eq. (48), we are keeping the corrections that come 
from |cs| ~ 1 to be more accurate. The region in the (|cs|,cii) plane which is consistent with 
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Figure 11: Analysis of the near de Sitter model with < and di < AHMpJ{M^H). We show 
the region of parameter space consistent with the WMAP upper limit /^/^'^' ^ 435 (95% CL), and 
the parameter constraints that must be satisfied in this model (eq. (95)). Notice that /^^''' must be 
positive in this model. Using loose constraints (taking K = 1 as defined below eq. (95)), a very small 
region of parameter space is consistent with WMAP. Using slightly more conservative constraints 
{K = 2), the entire parameter space is ruled out at 95% CL. 



the WMAP constraint f^f' < 435 (95% CL) from sec. 4.2 is shown in Fig. 11. Notice that 
in this case /^^'^' has to be positive. 

From eq. (44), the absolute value of the speed of sound is bounded to be 

Ic^l < I ■ (94) 

The consistency of the calculation (see [19]) forces the model parameters to satisfy the fol- 
lowing constraints: 

-di 1 , / 25 , , , ,2,4 rff \ 



4(^2 + 4)^/^ 2 V4608 



H H 



(1-6 






-di){d2 + dsy/' 


8|c, 


2 



> 1 . (95) 



The first constraint follows from the expression for the power spectrum in eq. (49) and the 
requirement that di < —4(^2 + d^Y^'^ in this model, which just comes from imposing that 
the dispersion relation at horizon crossing is dominated by the term is Csk. The second 
requirement comes from imposing that the cutoff of the model is much larger than Hubble. 

To make the constraints precise, let us introduce a "threshhold" parameter K > 1, and 
define the symbol ^ in eq. (95) to mean that the ratio between the left-hand and right-hand 
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sides must be > K. As shown in Fig. 11, if we take K — 1 (corresponding to a loose interpre- 
tation of the constraints) then a small region in parameter space is consistent with WMAP, 
whereas if we take K — 2 (corresponding to a slightly more conservative interpretation) then 
the entire parameter space is ruled out. 

As in the former cases, we have neglected metric fluctuations. Following [20, 10] it is 
straightforward to see that this is a good approximation for di < (Mpi/(30M))^. When this 
inequality is violated, as shown in [20, 10] the squared speed of sound of the fluctuations 
is still negative. In this case the model is still very similar to the one just shown, and, as we 
are seeing, it is generically very disfavored by the data. 

4.7 Near-de-Sitter models with < and \di\ < 4HM^J{M^H) 
In this last case, the bispectrum is given by: 



(/Z"') \ _ ( 1-040 0.9878 
{fNL°^') ) \ 0-1079 0.1613 



/ (5.179 X 10-) (™|a*D) 



by combining Eqs. (50) and (65). Notice that in this model H > and that, contrary to 
what we did in eq. (48), we are keeping the corrections that come from |cs| ~ 1 to be more 
accurate. 

Because the bispectrum depends on two combinations of parameters, we define a statis- 
tic which uses both /^^''' and f^^^^', by using eq. (89) with expectation values given by 
eq. (96). In Fig. 12, we show the region of parameter space which is consistent with WMAP, 
using this statistic. 

We have the following parameter constraints in this model: 

211/4 /^^-- . 



™ll(l + 
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log — A^|c,|(l + |c^|)— /I »1 



{d2 + d3y/m 2 ^ \^ 9 ^1 ' I 
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log(lf A^|c,|(l + |c2 




The first constraint follows from the requirement that < {4H Mp^) / {M^ H) in this model, 
together with the expression for the sound speed in eq. (47). The second constraint comes 
from imposing that the dispersion relation at horizon crossing is dominated by the term is Csk, 
while the third comes from imposing the cutoff scale to be larger than H. As in the previous 
subsection, we define the symbol ^ in eq. (97) by introducing a threshhold parameter K: we 
say that a; <^ y if the ratio {y/x) is > K. 
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Figure 12: Analysis of the near de Sitter model with < and \di\ < 4:HMp^/{M^H). Notice that 
this model has H > 0. We show the regions consistent with WMAP, defined using a statistic 
which includes /^^'^' and f'^^]^°^', and the regions that are disallowed by the constraints in eq. (97). 
Using loose constraints (taking K = 1 as defined below eq. (97)), there is a large region of parameter 
space which is consistent with WMAP, including the "branch" in the upper right part of the plot 
where f^]/^"^' is generated. Using somewhat more conservative constraints {K = 4), the entire 
parameter space is ruled out. 

If we use loose constraints {K = 1), then there is a large region of parameter space which 
is consistent with WMAP, including an extended "branch" which points toward the upper 
right corner in Fig. 12. This branch corresponds to models in which the contributions to the 
bispectrum from the operators 7r(97r)^ and (9^7r)(97r)^ have opposite signs and nearly cancel, 
so that the orthogonal shape is generated instead of the equilateral shape. Because our 
analysis includes f^^"^' , the branch is cut off and arbitrarily large values of the parameters 
are not allowed. 

If we use somewhat more conservative constraints {K = 4), then the entire parameter 
space is ruled out, including the branch. 

As in the former cases, we have neglected metric fluctuations. Following [20, 10] it is 
straightforward to see that this is a good approximation for + ^3 ^ (Mpi/(30M))^°. When 
this inequality is violated, as shown in [20, 10] the squared speed of sound cj. of the fluctuations 
is still negative. In this case the model is still very similar to the one just shown, and, as we 
are seeing, it is generically very disfavored by the data. 

5 Summary 

The non-Gaussianities produced by the most general single-fleld inflation with an approximate 
shift symmetry protecting the Goldstone boson are described by two independent parameters. 
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The resulting shape of the signal in Fourier space can vary from being peaked on equilateral 
triangles to being peaked on fiat triangles (where the two shortest sides are equal to half the 
longest one) . Two independent parameters control this bi-dimensional space of signatures and 
they both need to be constrained. We have shown that former analyses have been partially 
blind to a fraction of this parameter space. 

We have then shown that an analysis covering the whole of the parameter space can 
be performed in practice using two approximate factorizable templates. The first template 
is the standard equilateral one, which constrains the amount of signal that is peaked on 
equilateral configurations in Fourier space. We call the second template 'orthogonal' and we 
choose it to constrain the amount of signal that in Fourier space is peaked both on equilateral 
configurations and on flat triangular configurations (with opposite signs). This new shape is 
approximately orthogonal to the equilateral one using a 3D scalar product. 

We apply the optimal estimator for the two Jnl parameters that characterize the ampli- 
tude of the two independent templates to the WMAP 5-year data. We find no evidence of 
non-Gaussianities, and we obtain the following optimal constraints: 

-125 < ft^f- < 435 at 95% CL , (98) 
-369 < f°^^°^- < 71 at 95% CL . (99) 

The Corr(/^^''", Z^^**"^ ) in the Monte Carlo simulations is equal to 0.32. The constraint on 
/nl^ is milder than the one presented by the WMAP collaboration for the same set of data 
[22] . We beheve that previous analyses have underestimated the error bar on f^l'^ and that 
our constraints are optimal. 

Thanks to the effective field theory of inflation [3], we are able to map the constraints on 
the two Jnl parameters above into constraints on the coefficient of the interaction Lagrangian 
of the Goldstone boson. Under the assumption that the primordial density perturbations are 
generated by a single-field inflationary model where there is an approximate shift symmetry 
for the Goldstone boson, this mapping is unique and constrains all possible inflationary models 
of this kind. We stress again that we do not assume that the background solution is given 
by a fundamental scalar held: the Lagrangian for the fluctuations in terms of the Goldstone 
boson is independent of the details through which the background solution is generated. We 
are only assuming that there is one light degree of freedom playing a relevant role during the 
inflationary phase. We draw contour plots for the parameters of the Goldstone Lagrangian 
that arc constrained by our analysis. In particular, for one of them, the speed of sound of the 
inflaton fluctuations, we flnd that it has to be larger than: 

Cs > 0.011 at 95% CL , (100) 

or smaller than 

Cs < 10-' (d2 + dsf' . (101) 

In this last case the higher-derivative kinetic term is important at horizon crossing and the 
non-Gaussianities depend on other coefficients. In some region of the parameter space, con- 
sistent inflationary models have a negative squared speed of sound for the fluctuations at 
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horizon crossing. This leads to an exponential growth of the perturbations before horizon 
crossing and an associated increase in the level of non-Gaussianities. The analysis of the 
WMAP data shows that these models are practically ruled out at 95% CL. 

These error bars are expected to decrease by a factor of six in the next few years when 
the data from the Planck satellite will be available [49] . 
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A A different but equivalent definition of the estimator 

In the main text, we related /^^'^' and ^'^ and C3 in eq. (60) by noticing that the 

expectation value of our estimator was given by the 3D Fisher matrix between the shape of 
non-Gaussianity in single field infiation and the templates. If instead one wanted to define 
/^^''"(cs, 63) and (cs, C3) so that our template shapes best approximate the exact ex- 

pression for given values of Lagrangian parameters and 63, then one would instead use the 
following expression 



fNL i'^s, C3) 
'Orth 
NL 



JNL lCs,C3 



(102) 



-^cquil.,1 ■ -^equil. ,1 -^equil.,1 ' -^orthog.,1 \ ( -^cquil.,1 ' F 

-^equil. ,1 ' -^orthog. ,1 -^orthog.,1 ' -^orthog.,1 / \ -^orthog.,1 ' ^ 

1.056 1.283 \ / /Jf^")'(c,) 
-0.05721 -0.2663 J [ r^^ics.h) 



where the subscript 1 means that the shape is evaluated with /jvl = 1. This expression is 
in fact obtained by minimizing the Fisher Matrix distance defined as the norm of the shape 
F ~ -^template accordiug to the 3D scalar product. 



^^Here, for simplicity, we will refer only to the case where inflation happens in the way that in the text we 
are referring to as 'not-so-near-de-Sitter'. The discussion applies with trivial modifications to when inflation 
happens near-de-Sitter. 
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The reason why in the main text we use instead the definition in (60) is that, for what 
the data analysis is concerned, we are interested in the statistics, and this depends on the 
definition of the estimator. In fact, expression (60) represents the expectation value of the 
estimator in the presence of a given /^^^''^^ and /^^, in the approximation that the 3D scalar 
product is a good approximation to the 2D one. This is the correct relationship according to 
the definition of the estimator we use and that we explain in the main text, where we analyze 
the data for the equilateral and the orthogonal shape assuming, in each case, that the other 
shape gives zero contribution. 

However, this is not the only definition of the estimator one could have used. One could 
have in fact imagined to build an estimator for each shape that has zero response to the other 
shape. This different definition would give different results for the constraints on f^l'^' and 
Inl^^^ i b^^' goiiig to show, it would not change the statistics from which 

we deduce the constraints on the parameters of the Lagrangian. Let us see in detail how this 
happens, and to this purpose let us define the matrix F of the scalar products between the 
templates as 

pi j (-^cquil.,1 ■ -^equil.,l) (-fequil.,1 ' -^orthog.,l) j (103) 

\ (-^equil. ,1 ■ -^orthog.,l) (-^orthog. ,1 ' -^orthog.,l) / 

and the matrix M of the scalar products between the templates and the single field shapes as 

^ f (-^cquil.,1 ■ -F7r(9i7r)2) (-^cquil. , 1 ' -^tt^ ) \ ^-j^q^-j 
\ (-^orhtog., tt)'^) (-^orthog.,1 " -^7T^ ) / 

In the two equations above, as well as in the rest of this Appendix, the dot ■ is now meant 
to represent a 2D scalar product as defined in eq. (55). For simplicity, we call the induced 
byspectrum Bi by a shape Fi simply as Fj. 

Let us also define the vector of the two estimators 

(?equil. \ 
Jnl / 

where the / represents the fact that we have removed the normalization factors iVgquii. = Fn 
and Ai'orthog. = -^22 from eq. (74) and (75) respectively. We also define the vector / of the 
/atl's as 



/ = '^"3 . (106) 




The way we define the estimator E^^^ that we use in in our analysis is by normalizing the 
estimator E in the following way: 

1 Jcquil. 

^(1) = I W^JnL I .^Q7^ 
I J_ /-orthog. ' ' ^ ^ 

F22 ^ 
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Since at the moment we are not seeing any evidence for a non-zero value for any of the /nlS, 
the estimator E'^^^ analyzes the data for each shape assuming that the other shape is not 
present. The covariance C^^^ for this estimator is given by the matrix 

= ^if- , (108) 

V F11F22 F22 ) 

In the case of the estimator it is straightforward to see that the ^ statistics is given by 

X(l)(/) = (/ - fwMApfM^F-'M{f - /WMAP) , (109) 

where /wmap is the result of / deduced from the WMAP data, and we have made use of 
the matrix M to translate the values of f^l^^' and f^^i"^' into values for fP^'^^ and f^\. 
The matrix M we use in the analysis is analogous to the one defined in eq. (104), with the 
2D scalar product replaced with 3D scalar products. We believe that this is a reasonable 
approximation, at order 10% level, given that the templates we use are very similar to the 
actual shape. This level of accuracy is in particular justified by the current obervational 
status, where we do not have a detection, but only constraints. 

Now, we could instead have defined for our analysis the estimator in an alternative way, 
that we call E^'^^: 

^(2) ^ _ (^^0) 

This definition is also commonly used in the hterature. With this definition, there is no cross- 
response between the equilateral and the orthogonal shape, and therefore the expectation 
value of the estimator in the presence of an hypothetical non-Gaussianity in the sky induced 
by single field infiation with a given Cg and £3, is given by eq. (102) with the 3D scalar products 
replaced with 2D ones. In this case, the estimator covariance C^"^^ is simply given by 

C(2) = , (111) 

and the is given by 

X(2)(/) = (/ - /wMAp)^M^F-iF-iFM(/ - /wmap) = (112) 
(/ - fwMApfM^F-^M{f - /wmap) . 

where we have used that the matrix in eq. (102) is equal to p-^M and we have made the same 
approximation as before in approximating 2D scalar products with 3D scalar products. The 
defined for the two estimators E^^^ and E^'^\ using respectively the mappings of eq. (60) 
and eq. (102), are equivalent. This means that while the constraints we get on f^l'^' and 
fNL^°^ (io depend on the choice of the estimator (and we explain in the main part of the text 
why we prefer the first definition of the estimators) , the constraints on the parameters of the 
Lagrangian, do not depend on these two choices. A very similar argument to the one above 
shows that our constraints on single-field infiation would be unchanged if we used a different 
value of c to define /nl""^' in eq. (58). 
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B A more accurate template 



We have found that a different choice of templates with respect to the one we make in the 
main text allows for a better approximation of the bi-dimensional space of non-Gaussianities 
produced in single field inflation. This is achieved if one substitutes the orthogonal template 
we have used in the main part of the paper with the following one, which is again similar to 
the orthogonal shape (see upper- right panel of Fig. 1): 



-^orthog. (2) — /atl""^ ■ 

Here p ~ 8.52 , and 



A(/ci,/c2,/c3) r{ki,k2,h) 

i^+P) ,.3 ,.3 ,.3 P- 



12 3 '^l''^2 3 



(113) 



A(A;i, k2, ks) = {h - ki){kt - k2){kt - k^) , (114) 
with kt — ki + k2 + ks, and 

r(A;i, k2, ks) = ^(fcifc2 + k2k3 + k^k,) - ^{kl + ki + A;|) . (115) 

The value of p is chosen in such a way that the template Forthog. (2) has the useful property of 
being orthogonal to the template -Fequu, in the sense of the 3D scalar product between shapes 
of eq. (23). In practice, Fcquu is a good template for the equilateral shape with 63 ~ 0, while 
Forthog. is a good template for the orthogonal shape £3 ~ —5.4. 

Similarly to what we did in the main part of the text, we define the template for single-field 
inflation to be: 

-^singl. field(2)(A;i, A;2, /C3, C^, C3) = /^T^'(Cs, C3)Fequii.(/Cl, k2, k3)+f^^''^'^'^\cs, C3)Forthog. {2)(ki, h, h) 

(116) 

The analogous of eq. (60), now becomes 

' \ \ orthog. (2) --f orthog. (2) / \ -''orthog. (2) '-''orthog. (2) / / /nl-A| = 1 

(117) 



1.040 1.210 \ / (c) 
-0.03951 -0.1757 I /^Kccg) 



where f^^''^'"^ (c^) and f^l'""^ (c^jCs) are given in eq. (21). 

Quite remarkably, the scalar product of this template with the exact shape (in the case 
not-near-de-Sitter) for any value of Cg and £3 is always larger than 0.99 (the minimum is at 
C3 ~ —5.4). This tells us that these two independent templates are effectively covering all 
the parameter space of shapes generated in single-field inflation. The reason why this second 
orthogonal template does such a better job than the one we use in the text is that in the 
squeezed limit k\ <S k2,ks, Forthog. (2)(^i, ^2, ^3) ~ 0{l/ki) as the exact shape, while the 
template we use in the text goes as 0{l/kl). The amount of signal contained in this limit 
is not very large, and this is why the orthogonal template we use in the main text gives still 
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a good approximation to the exact shape. We have found that implementing this template 
numerically is quite difficult, and therefore, given than the orthogonal template we present 
in the main text does not create such numerical difficulties and it is still accurate enough, we 
have decided to perform the analysis of the data with the Forthog. presented in the text. 
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